3.136 \(\int (g+h x)^m (a+c x^2)^p (d+e x+f x^2) \, dx\)

Optimal. Leaf size=420 \[ -\frac{\left (a+c x^2\right )^p (g+h x)^{m+1} \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{\frac{\sqrt{-a} h}{\sqrt{c}}+g}\right )^{-p} F_1\left (m+1;-p,-p;m+2;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right ) \left (a f h^2 (m+1)-c \left (2 f g^2 (p+1)-h (m+2 p+3) (e g-d h)\right )\right )}{c h^3 (m+1) (m+2 p+3)}-\frac{\left (a+c x^2\right )^p (g+h x)^{m+2} \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{\frac{\sqrt{-a} h}{\sqrt{c}}+g}\right )^{-p} (2 f g (p+1)-e h (m+2 p+3)) F_1\left (m+2;-p,-p;m+3;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )}{h^3 (m+2) (m+2 p+3)}+\frac{f \left (a+c x^2\right )^{p+1} (g+h x)^{m+1}}{c h (m+2 p+3)} \]

[Out]

(f*(g + h*x)^(1 + m)*(a + c*x^2)^(1 + p))/(c*h*(3 + m + 2*p)) - ((a*f*h^2*(1 + m) - c*(2*f*g^2*(1 + p) - h*(e*
g - d*h)*(3 + m + 2*p)))*(g + h*x)^(1 + m)*(a + c*x^2)^p*AppellF1[1 + m, -p, -p, 2 + m, (g + h*x)/(g - (Sqrt[-
a]*h)/Sqrt[c]), (g + h*x)/(g + (Sqrt[-a]*h)/Sqrt[c])])/(c*h^3*(1 + m)*(3 + m + 2*p)*(1 - (g + h*x)/(g - (Sqrt[
-a]*h)/Sqrt[c]))^p*(1 - (g + h*x)/(g + (Sqrt[-a]*h)/Sqrt[c]))^p) - ((2*f*g*(1 + p) - e*h*(3 + m + 2*p))*(g + h
*x)^(2 + m)*(a + c*x^2)^p*AppellF1[2 + m, -p, -p, 3 + m, (g + h*x)/(g - (Sqrt[-a]*h)/Sqrt[c]), (g + h*x)/(g +
(Sqrt[-a]*h)/Sqrt[c])])/(h^3*(2 + m)*(3 + m + 2*p)*(1 - (g + h*x)/(g - (Sqrt[-a]*h)/Sqrt[c]))^p*(1 - (g + h*x)
/(g + (Sqrt[-a]*h)/Sqrt[c]))^p)

________________________________________________________________________________________

Rubi [A]  time = 0.604235, antiderivative size = 417, normalized size of antiderivative = 0.99, number of steps used = 6, number of rules used = 4, integrand size = 27, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.148, Rules used = {1654, 844, 760, 133} \[ -\frac{\left (a+c x^2\right )^p (g+h x)^{m+1} \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{\frac{\sqrt{-a} h}{\sqrt{c}}+g}\right )^{-p} F_1\left (m+1;-p,-p;m+2;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right ) \left (a f h^2 (m+1)+c h (m+2 p+3) (e g-d h)-2 c f g^2 (p+1)\right )}{c h^3 (m+1) (m+2 p+3)}-\frac{\left (a+c x^2\right )^p (g+h x)^{m+2} \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{\frac{\sqrt{-a} h}{\sqrt{c}}+g}\right )^{-p} (2 f g (p+1)-e h (m+2 p+3)) F_1\left (m+2;-p,-p;m+3;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )}{h^3 (m+2) (m+2 p+3)}+\frac{f \left (a+c x^2\right )^{p+1} (g+h x)^{m+1}}{c h (m+2 p+3)} \]

Antiderivative was successfully verified.

[In]

Int[(g + h*x)^m*(a + c*x^2)^p*(d + e*x + f*x^2),x]

[Out]

(f*(g + h*x)^(1 + m)*(a + c*x^2)^(1 + p))/(c*h*(3 + m + 2*p)) - ((a*f*h^2*(1 + m) - 2*c*f*g^2*(1 + p) + c*h*(e
*g - d*h)*(3 + m + 2*p))*(g + h*x)^(1 + m)*(a + c*x^2)^p*AppellF1[1 + m, -p, -p, 2 + m, (g + h*x)/(g - (Sqrt[-
a]*h)/Sqrt[c]), (g + h*x)/(g + (Sqrt[-a]*h)/Sqrt[c])])/(c*h^3*(1 + m)*(3 + m + 2*p)*(1 - (g + h*x)/(g - (Sqrt[
-a]*h)/Sqrt[c]))^p*(1 - (g + h*x)/(g + (Sqrt[-a]*h)/Sqrt[c]))^p) - ((2*f*g*(1 + p) - e*h*(3 + m + 2*p))*(g + h
*x)^(2 + m)*(a + c*x^2)^p*AppellF1[2 + m, -p, -p, 3 + m, (g + h*x)/(g - (Sqrt[-a]*h)/Sqrt[c]), (g + h*x)/(g +
(Sqrt[-a]*h)/Sqrt[c])])/(h^3*(2 + m)*(3 + m + 2*p)*(1 - (g + h*x)/(g - (Sqrt[-a]*h)/Sqrt[c]))^p*(1 - (g + h*x)
/(g + (Sqrt[-a]*h)/Sqrt[c]))^p)

Rule 1654

Int[(Pq_)*((d_) + (e_.)*(x_))^(m_.)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> With[{q = Expon[Pq, x], f = Coeff
[Pq, x, Expon[Pq, x]]}, Simp[(f*(d + e*x)^(m + q - 1)*(a + c*x^2)^(p + 1))/(c*e^(q - 1)*(m + q + 2*p + 1)), x]
 + Dist[1/(c*e^q*(m + q + 2*p + 1)), Int[(d + e*x)^m*(a + c*x^2)^p*ExpandToSum[c*e^q*(m + q + 2*p + 1)*Pq - c*
f*(m + q + 2*p + 1)*(d + e*x)^q - f*(d + e*x)^(q - 2)*(a*e^2*(m + q - 1) - c*d^2*(m + q + 2*p + 1) - 2*c*d*e*(
m + q + p)*x), x], x], x] /; GtQ[q, 1] && NeQ[m + q + 2*p + 1, 0]] /; FreeQ[{a, c, d, e, m, p}, x] && PolyQ[Pq
, x] && NeQ[c*d^2 + a*e^2, 0] &&  !(EqQ[d, 0] && True) &&  !(IGtQ[m, 0] && RationalQ[a, c, d, e] && (IntegerQ[
p] || ILtQ[p + 1/2, 0]))

Rule 844

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Dist[g/e, Int[(d
+ e*x)^(m + 1)*(a + c*x^2)^p, x], x] + Dist[(e*f - d*g)/e, Int[(d + e*x)^m*(a + c*x^2)^p, x], x] /; FreeQ[{a,
c, d, e, f, g, m, p}, x] && NeQ[c*d^2 + a*e^2, 0] &&  !IGtQ[m, 0]

Rule 760

Int[((d_) + (e_.)*(x_))^(m_)*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> With[{q = Rt[-(a*c), 2]}, Dist[(a + c*x^
2)^p/(e*(1 - (d + e*x)/(d + (e*q)/c))^p*(1 - (d + e*x)/(d - (e*q)/c))^p), Subst[Int[x^m*Simp[1 - x/(d + (e*q)/
c), x]^p*Simp[1 - x/(d - (e*q)/c), x]^p, x], x, d + e*x], x]] /; FreeQ[{a, c, d, e, m, p}, x] && NeQ[c*d^2 + a
*e^2, 0] &&  !IntegerQ[p]

Rule 133

Int[((b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_)*((e_) + (f_.)*(x_))^(p_), x_Symbol] :> Simp[(c^n*e^p*(b*x)^(m +
 1)*AppellF1[m + 1, -n, -p, m + 2, -((d*x)/c), -((f*x)/e)])/(b*(m + 1)), x] /; FreeQ[{b, c, d, e, f, m, n, p},
 x] &&  !IntegerQ[m] &&  !IntegerQ[n] && GtQ[c, 0] && (IntegerQ[p] || GtQ[e, 0])

Rubi steps

\begin{align*} \int (g+h x)^m \left (a+c x^2\right )^p \left (d+e x+f x^2\right ) \, dx &=\frac{f (g+h x)^{1+m} \left (a+c x^2\right )^{1+p}}{c h (3+m+2 p)}+\frac{\int (g+h x)^m \left (-h^2 (a f (1+m)-c d (3+m+2 p))-c h (2 f g (1+p)-e h (3+m+2 p)) x\right ) \left (a+c x^2\right )^p \, dx}{c h^2 (3+m+2 p)}\\ &=\frac{f (g+h x)^{1+m} \left (a+c x^2\right )^{1+p}}{c h (3+m+2 p)}+\frac{\left (e h-\frac{2 f g (1+p)}{3+m+2 p}\right ) \int (g+h x)^{1+m} \left (a+c x^2\right )^p \, dx}{h^2}-\frac{\left (a f h^2 (1+m)-2 c f g^2 (1+p)+c h (e g-d h) (3+m+2 p)\right ) \int (g+h x)^m \left (a+c x^2\right )^p \, dx}{c h^2 (3+m+2 p)}\\ &=\frac{f (g+h x)^{1+m} \left (a+c x^2\right )^{1+p}}{c h (3+m+2 p)}+\frac{\left (\left (e h-\frac{2 f g (1+p)}{3+m+2 p}\right ) \left (a+c x^2\right )^p \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p}\right ) \operatorname{Subst}\left (\int x^{1+m} \left (1-\frac{x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^p \left (1-\frac{x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^p \, dx,x,g+h x\right )}{h^3}-\frac{\left (\left (a f h^2 (1+m)-2 c f g^2 (1+p)+c h (e g-d h) (3+m+2 p)\right ) \left (a+c x^2\right )^p \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p}\right ) \operatorname{Subst}\left (\int x^m \left (1-\frac{x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^p \left (1-\frac{x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^p \, dx,x,g+h x\right )}{c h^3 (3+m+2 p)}\\ &=\frac{f (g+h x)^{1+m} \left (a+c x^2\right )^{1+p}}{c h (3+m+2 p)}-\frac{\left (a f h^2 (1+m)-2 c f g^2 (1+p)+c h (e g-d h) (3+m+2 p)\right ) (g+h x)^{1+m} \left (a+c x^2\right )^p \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} F_1\left (1+m;-p,-p;2+m;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )}{c h^3 (1+m) (3+m+2 p)}+\frac{\left (e h-\frac{2 f g (1+p)}{3+m+2 p}\right ) (g+h x)^{2+m} \left (a+c x^2\right )^p \left (1-\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} \left (1-\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )^{-p} F_1\left (2+m;-p,-p;3+m;\frac{g+h x}{g-\frac{\sqrt{-a} h}{\sqrt{c}}},\frac{g+h x}{g+\frac{\sqrt{-a} h}{\sqrt{c}}}\right )}{h^3 (2+m)}\\ \end{align*}

Mathematica [F]  time = 1.19881, size = 0, normalized size = 0. \[ \int (g+h x)^m \left (a+c x^2\right )^p \left (d+e x+f x^2\right ) \, dx \]

Verification is Not applicable to the result.

[In]

Integrate[(g + h*x)^m*(a + c*x^2)^p*(d + e*x + f*x^2),x]

[Out]

Integrate[(g + h*x)^m*(a + c*x^2)^p*(d + e*x + f*x^2), x]

________________________________________________________________________________________

Maple [F]  time = 0.796, size = 0, normalized size = 0. \begin{align*} \int \left ( hx+g \right ) ^{m} \left ( c{x}^{2}+a \right ) ^{p} \left ( f{x}^{2}+ex+d \right ) \, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((h*x+g)^m*(c*x^2+a)^p*(f*x^2+e*x+d),x)

[Out]

int((h*x+g)^m*(c*x^2+a)^p*(f*x^2+e*x+d),x)

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (f x^{2} + e x + d\right )}{\left (c x^{2} + a\right )}^{p}{\left (h x + g\right )}^{m}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((h*x+g)^m*(c*x^2+a)^p*(f*x^2+e*x+d),x, algorithm="maxima")

[Out]

integrate((f*x^2 + e*x + d)*(c*x^2 + a)^p*(h*x + g)^m, x)

________________________________________________________________________________________

Fricas [F]  time = 0., size = 0, normalized size = 0. \begin{align*}{\rm integral}\left ({\left (f x^{2} + e x + d\right )}{\left (c x^{2} + a\right )}^{p}{\left (h x + g\right )}^{m}, x\right ) \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((h*x+g)^m*(c*x^2+a)^p*(f*x^2+e*x+d),x, algorithm="fricas")

[Out]

integral((f*x^2 + e*x + d)*(c*x^2 + a)^p*(h*x + g)^m, x)

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((h*x+g)**m*(c*x**2+a)**p*(f*x**2+e*x+d),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (f x^{2} + e x + d\right )}{\left (c x^{2} + a\right )}^{p}{\left (h x + g\right )}^{m}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((h*x+g)^m*(c*x^2+a)^p*(f*x^2+e*x+d),x, algorithm="giac")

[Out]

integrate((f*x^2 + e*x + d)*(c*x^2 + a)^p*(h*x + g)^m, x)